Photoassociative Production and Trapping of Ultracold KRb Molecules 
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We have produced ultracold heteronuclear KRb molecules by the process of photoassociation in 
a two-species magneto-optical trap. Following decay of the photoassociated KRb*, the molecules 
are detected using two-photon ionization and time-of- flight mass spectroscopy of KRb + . A portion 
of the metastable triplet molecules thus formed are magnetically trapped. Photoassociative spectra 
down to 9f cm -1 below the K(4s) + Rb (5pi/2) asymptote have been obtained. We have made 
assignments to all eight of the attractive Hund's case (c) KRb* potential curves in this spectral 
region. 
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Following in the footsteps of ultracold atoms, the study 
of ultracold molecules has seen significant progress in re- 
cent years Homonuclear molecules have been pro- 
duced at ultracold temperatures, here defined as T < 1 
mK, by photoassociating ultracold atoms |3j. In a 
few cases, the resulting molecules have been subsequently 
trapped |4J,|5|. Feshbach resonances in even colder atomic 
samples have yielded very weakly bound homonuclear 
molecules 0, in some cases under quantum degenerate 
conditions. Heteronuclear polar molecules have attracted 
particular interest because their permanent dipole mo- 
ments can be controlled with external fields. At ultra- 
cold temperatures, the resulting long-range, anisotropic 
dipole-dipole interactions are expected to lead to new 
phenomena, such as anisotropic collisions Q and novel 
quantum degenerate behavior Jj|, 0, [^j . Heteronuclear 
photoassociative spectroscopy [ll| is also of interest, and 
applications of ultracold polar molecules to quantum 
computing [l^ and fundamental symmetries [l3j are also 
being considered. 

In this Letter, we report on the production, detection, 
magnetic trapping, and photoassociation spectroscopy of 
KRb molecules. This work builds on a variety of other 
recent advances, particularly the recent production and 
spectroscopic study of ultracold polar RbCs molecules by 
Kerman and coworkers |l4i fl5| , using photoassociation 
(PA) in overlapping Rb and Cs magneto-optical traps 
(MOTs). Similar PA spectroscopy has been performed 
with the heteronuclear (but nonpolar) 6 Li 7 Li molecule 
(lif. In other related work, observations of ultracold 
NaCs+ [13, RbCs H3, and KRb [HI have been reported 
in two-species MOTs, although the formation mecha- 
nisms have not been identified. Heteronuclear Feshbach 
resonances have been recently reported in Li+Na [2fJ and 
K+Rb [2i|, although the resulting molecules were not 
directly observed. At temperatures above 10 mK, po- 



lar molecules have been successfully cooled and trapped. 
CaH has been buffer-gas cooled and magnetically trapped 
[2^|. while ND3 has been electrostatically slowed and 
loaded into electric traps [23H24| . 

Photoassociation, as shown in Figure 1, involves the 
free-bound absorption of a photon by a colliding pair of 
atoms Q • It is a particularly powerful spectroscopic 
technique at /iK temperatures, where only a few partial 
waves contribute, as the centrifugal barriers of collision 
are significantly greater than the thermal kinetic energy 
for quite low collisional angular momenta (£ < 3). Thus 
the bound states formed are also limited to low rotational 
angular momenta ( J < 4 ). 

In our experiments, electronically excited KRb* 
molecules are produced by PA in overlapping K and Rb 
MOTs, followed by radiative stabilization via decay to 
the ground X 1 T J + state or the metastable et 3 £ + state. 
Dark-spot MOTs are used to achieve high atomic densi- 
ties, and ionization detection is employed for improved 
sensitivity. We are able to observe trapping of ultracold 
heteronuclear molecules for the first time, by confining 
triplet KRb molecules in the inhomogeneous magnetic 
field of the MOT. We also present high-resolution pho- 
toassociative spectra down to 91 cm" 1 below the K(4s) 
+ Rb (5p!/ 2 ) asymptote. In these spectra we have iden- 
tified transitions to all eight of the attractive long-range 
potentials converging to 4s + bpi/2 and 4s + 5p 3 / 2 . 

Our 39 K/ 85 Rb dual species MOT uses the same 
apparatus previously used for trapping 39 K [2(| with the 
addition of trap and repump beams for 85 Rb. Frequently 
the K and Rb diode laser beams can share optics, because 
of the closeness of the resonance energies (see Fig. 1). 
Indeed it is this near degeneracy which makes the KRb 
Ce values especially large ^jj . Normally we operate each 
of the atom traps as a "dark-spot" MOT j^lj , with typical 
densities of 10 n /cm 3 for Rb and 3 x 10 10 /cm 3 for K. The 
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FIG. 1: The photo-association process (PA) and possible sin- 
glet (X J E + ) and triplet (a E + ) molecule formation processes 
in KRb based on the ab initio potential energy curves of Ref. 
[29|| . In (b), the long-range states converging to 4s + 5p./ arc 
shown; the state designations consist of an ordinal number- 
ing followed by the Hund's case (c) symmetry in parenthe- 
ses, expressed as Q ± . The long-range states converging to 
4pj + 5s 1 / 2 (not shown) are all repulsive. 

two MOTs are carefully aligned to maximize overlap of 
the two atomic clouds. 

To achieve sensitive and species-selective detection we 
employ two-photon ionization at 602.5 nm followed by 
time-of-flight mass spectroscopy, an example of which is 
shown in the inset of Fig. 2. The ionizing laser pulse (1.5 
mJ, ^10 ns) is focused to provide a typical peak intensity 
of 7 x 10 6 W/cm 2 . These pulses are generated using a 
doubled Nd:YAG laser at 532nm to pump a pulsed dye 
laser at 10Hz. Ion signals can also be obtained using a 532 
nm laser directly. A Channcltron is used to detect atomic 
and molecular ions. We obtain PA spectra by tuning a 
cw tunable PA laser (Coherent 829-29, typically > 400 
mW) while detecting the KRb + time-of-flight peak. 

We obtain typical signals of 10-60 KRb + ions per 
laser shot when the PA laser is tuned to a strong KRb 
resonance. Surprisingly, these signals are comparable 
to the RbJ ion signals observed when photoassociat- 
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FIG. 2: Main panel: Time evolution of the d9 K 85 Rb + signal 
as a function of the delay of the ionization pulse relative to 
the turnoff of the PA laser, with the PA laser set to a 3(0~) 
state resonance as indicated in Fig. EJa). PA laser detuning 
Apa is measured relative to K(4s, F=l) + Rb(5pi/ 2 , F=2). 
Inset: Typical time-of-flight mass spectrum, with the PA laser 
tuned to the same resonance. 



ing on a strong Rb2 line. The number of colliding 
atom pairs available at the Condon radius, where the 
PA laser matches the ground-excited free-bound transi- 
tion, strongly favors Rb2. Its resonant-dipole 1/i? 3 po- 
tentials permit photoassociation at much longer range 
(e.g., 100 ao) compared to KRb where the van der Waals 
1/i? 6 potential results in shorter-range excitation (e.g., 
30 ao). However, the Franck-Condon factors for decay 
to the 1/i? 6 ground-state potentials strongly favor KRb. 
Therefore we expect to produce less KRb* in the PA step, 
but with more efficient decay to a bound state. This is 
reflected in the barely detectable KRb trap loss signal 
we observe in contrast to losses of 10-50% for Rb2. Also, 
the relaxed heteronuclear selection rules allow more PA 
pathways in KRb compared to Rb2 ^3 • The signal sizes 
may also be influenced by different detection efficiencies 
for the KRb and Rb2. Finally, the signal sizes may be en- 
hanced by magnetic trapping, discussed below. Ignoring 
this trapping and assuming 100% ionization efficiency for 
KRb and 50% ion detection efficiency, our maximum sig- 
nal size would correspond to a production rate of about 
4 x 10 4 KRb molecules/s. 

We observe trapping of triplet KRb molecules in the 
magnetic field of the MOT coils |5J. These molecules 
must be in the triplet a state since the singlet molecules 
have zero magnetic moment. We can investigate the trap- 
ping by turning off the PA laser while leaving the MOT 
lasers and magnetic field on, then firing the detection 
laser pulse after a variable delay time. Figure 2 shows a 
decay curve observed as a function of this delay, with the 
laser tuned to one of the stronger KRb resonances (the 
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FIG. 3: (a) A PA spectrum for 39 K 85 Rb showing eight vi- 
brational bands as our PA laser is scanned as indicated below 
the K(4s) + Rb(5pi/2) asymptote. The tentative state desig- 
nation is given for each band using the nomenclature of Fig 
1(b). The rotational spacings vary widely, with smaller spac- 
ings corresponding to longer-range states. This spectrum in- 
cludes transitions to seven of the eight states correlating with 
the J = 1/2 and J — 3/2 asymptotes. All these states can be 
reached by dipole-allowed transitions from colliding ground 
state atoms (Q — ± , 1). (b) A high-resolution PA spectrum 
showing the rotational structure as the PA laser detuning is 
scanned through a 2(0 + ) state resonance below the K(4s) + 
Rb(5pi/2) asymptote. The asterisks(*) indicate "hyperfine 
ghosts" due to a small population in the "bright" F=2 state 
of K(4s) in our dark-spot MOT. 



3(CP) transition marked with an arrow in Fig. 3(a)). The 
background of 1.7 ions/shot due to KRb formed in the 
absence of the PA laser has been subtracted from 
this decay curve. A significant fraction of the molecules 
survive long past the ballistic decay time, which is about 
3 ms for our 1 mm diameter detection beam, assuming a 
temperature of ^300 fiK. The trapped molecules appear 
to decay non-exponentially, with a small signal still ob- 
servable after 150 ms. Because the atoms in the MOT 
are unpolarized, we expect the molecules to have a wide 
distribution of projection quantum numbers (mj ormy, 
depending on the coupling scheme) , giving rise to a large 
spread in the effective trapping potential. For a molecule 
with the largest possible magnetic dipole moment projec- 
tion of two Bohr magnetons (2/j.b — 9e^/w e ), the hori- 
zontal trapping potential reaches a depth of 500 /iK with 
a radius of 2.5 mm. The vertical potential is slightly dis- 
torted by gravity, but the magnetic trapping force is 2.5 
times larger than the gravitational force. Molecules with 
smaller magnetic dipole moment projections will experi- 



ence weaker trapping, and can be expected to be con- 
fined in a larger trap volume, with a higher escape rate 
and lower detection efficiency. At least 2/3 of the total 
molecules formed will be untrapped or anti-trapped. 

Examples of our heteronuclear PA spectra are given 
in Figure 3. The observed linewidths vary from about 
35-500 MHz, primarily due to unresolved internal struc- 
ture. The excellent signal-to-noise ratio makes the spec- 
tra readily assignable. At these low J values, rotational 
series such as the one shown in Fig. 3(b) are fit precisely 
by a linear plot of energy versus (J(J+ 1) — O 2 ), yielding 
a rotational constant B v without observable centrifugal 
distortion. The four ± states (two of each), the three 

= 1 states, and the sole ft = 2 state are readily as- 
signed in this way. Our ± assignments are preliminary 
because of extensive perturbations |2^| . These effects are 
evident in Fig. 4, which plots the vibrational energies 
and B v values for one of the + states. We make nominal 
assignments by associating the state having the smaller 
Cg coefficient with the shorter outer turning point and 
therefore the larger rotational constant. Thus the se- 
ries in Figure 4 is assigned to the + state at the K(4s) 
+ Rb(5p!/ 2 ) asymptote, i.e. the 2(0 + ) state of Figure 
1(b). Figure 3(b) in particular corresponds to the ninth 
observed vibrational level in the Figure 4 series. There 
are eight vibrational levels in Figure 3a, assigned as or 

1 states depending on their lowest rotational quantum 
number ( J = or 1). One of the levels corresponds to 
the 11th level in the 2(0 + ) series in Fig. 0] Two other 
levels must correspond to 2(0~), with a vibrational spac- 
ing comparable to 2(0 + ) at this detuning. The final two 
levels, with larger vibrational spacings, are assigned 
to the 3(0~) and 3(0 + ) states, respectively, in order of 
increasing B v . 

In a similar way, among the three = 1 levels, one 
(with B v = 8.49 x 10~ 3 cm -1 ) fits well in a long vibra- 
tional series, similar to the 2(0 + ) series in Fig. 4, and is 
assigned to the 2(1) state. The other two levels cannot be 
in the 2(1) series and are assigned to the 3(1) and the 4(1) 
states based on their relative rotational constants (10.39 
and 10.93 x 10~ 3 cm -1 , respectively). Since four vibra- 
tional levels of the 1(2) state have also been observed at 
other detunings, we have been able to observe all eight 
of the Hund's case (c) electronic states correlating to the 
K(4s) + Rb(5p!/ 2 , 3/2) asymptotes, as shown in Fig 1(b). 
Further discussion of these assignments will be published 
elsewhere |28|. 

Due to our ion extraction configuration, the molecules 
are formed in a dc field of ~160V/cm. We observe no 
Stark effect at our spectral resolution of ^0.001 cm -1 , 
consistent with the small dipole moment (<0.04 eao) pre- 
dicted for KRb* E|j«J III |I2 This differs from the 
RbCs case of Ref. |l4j |. where both the field (390 V/cm) 
and the dipole moment (^0.5 eao) are larger. We note 
that the dipole moments of the AT 1 !]" 1 " and a 3 £ + states 
of KRb at their equilibrium separations are -0.30(2) and 
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FIG. 4: (a) The rotational constants B v in cm -1 and (b) the 
detuning (binding energy) to the 1/3 power for levels of the 
2(0 + ) states converging to the K(4s)+Rb(5pi/2) asymptote 
versus relative vibrational quantum number (absolute values 
unknown). For a pure RT 6 Ion g-range potential , th ese plots 
should be linear in the absence of perturbations |34l |35| . 



-0.02(1) eao, respectively [33j . 

Assuming our assignments are correct for the 2(0 + ), 
2(0~), and 2(1) states, one can estimate the correspond- 
ing effective C 6 values E3 as (1.06 ± 0.06) x 10 5 , 
(1.01 ± 0.04) x 10 5 and (1.01 ± 0.05) x 10 5 atomic units 
(a.u.), respectively. Note that Cq is a weak function of 
R between very long-range (spin-orbit 3> dispersion) and 
intermediate range (dispersion ^> spin-orbit) — the calcu- 
lations of [3||, for example, vary from 1.0 x 10 5 a.u. at 
very long range to 3.1 x 10 5 a.u. at intermediate range 
for the 2(0 + ) state. Model calculations to better define 
the dispersion interaction are underway. 

It is important to consider the effects of symmetry on 
molecule formation. Some relevant transition dipole mo- 
ment functions are available 0] . For levels near dissoci- 
ation, all of the excited states we produce decay with 
at least 78% probability to the a 3 S + state. In fact, 
the 2(0"), 3(0"), and 1(2) states radiate 100% to this 
triplet state. Thus, since we observe KRb + from 0" and 
2 states, we know our ion detection is sensitive to triplet 
molecules. However, since even the + and 1 states ra- 
diate primarily to triplet states at long range, we do not 
yet know if our ionization detection is sensitive to singlet 
molecules. We hope that in future experiments where 
both singlet and triplet molecules are produced but only 
triplets are trapped, we can determine the relative sensi- 
tivity for singlet and triplet molecule detection. 

In summary, we have significantly advanced the study 
of ultracold polar molecules, trapping them for the first 
time. KRb molecules are produced by high-resolution 
photoassociation and detected by ionization. We have ex- 
tensively explored the PA spectrum, identifying all eight 
of the electronic states expected near the K(4s) + Rb(5p) 
limits. 

We gratefully acknowledge support from the National 
Science Foundation and the University of Connecticut 
Research Foundation and helpful discussions with War- 
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